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Abstract 



We report the results of numerical calculations of rung-rung current cor- 
relations on a 2-leg t-J ladder with J/t = 0.35 for dopings x = 0.125 and 
x = 0.19. We find that the amplitude of these correlations decays exponen- 
tially. We argue that this can be understood within a bosonization framework 
in terms of the pinned phase variables associated with a C1SO phase with 
d x 2„j / 2-like power law pairing correlations. 
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Variational flux-phase states are known to provide interesting low energy estimates for 
the ground state of the 2D Heisenberg antiferromagnet. In addition, recent calculations 
have found signatures of a staggered-flux phase in current vorticity correlations in lightly- 
doped t-J clusters. In particular, Ivanov et. al. [|I| found evidence for a power law decay 
of the staggered current vorticity for a projected c?-wave variational wave function on a 
10 x 10 lattice and Leung [0 has reported staggered vorticity correlations from the exact 
diagonalization of a 32-site t-J cluster with two holes. Prompted by these results, we have 
carried out density matrix renormalization group (DMRG) calculations of the rung-rung 
current correlations on doped 2-leg t-J ladders. As is known, near half-filling an isotropic 
t-J ladder, with J/t in the physical range appropriate to the cuprates, retains a spin gap 
and exhibits power law <i x 2_j / 2-like pair field and 4fcp-CDW correlations 0. Thus, it is of 
interest to explore whether there may also be staggered current correlations on the 2-leg 
ladder. 

Indeed, as discussed below, the DMRG calculations show that there are oscillating rung- 
rung current correlations. However, we find that they are incommensurate at a finite doping, 
relatively weak, and decay exponentially on a scale of several lattice spacings. Complement- 
ing the numerical results, we discuss a bosonization treatment of the 2-leg ladder. Here, one 
finds that the rung current depends exponentially on various boson phase fields. We argue 
that in a C1SO phase which has a spin-gapped ground state with power law d x 2_ y 2-Y\ke pair- 
ing correlations, a phase variable appearing in the rung-current is dual to a pinned variable 
and that this leads to an exponential decay of the rung-rung current correlations, as seen 
numerically. 

Here we will consider an isotropic 2-leg t-J ladder described by the Hamiltonian 

H = 1 £ (4+lAaCiAa + h.C.) - tj^ {c\ 2a C ila + h.C.) 
i,X,a i,a 

— J ^2 ■ Si + w — lX * +1A j — J ( Sji ■ & 

Here Ci\ a destroys an electron on rung i and leg A = 1,2 with spin a =|, { and Si\ = 
c\ x a/2 Ci\. The Hilbert space excludes all states with doubly-occupied sites. The rung 
current operator for the j th rung can be written as 

= ti^2 ( C i2a C jla _ c ]la C j2a) = H ^ (^jea^joa ~ ^joa^jea) (2) 
a a 
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Here the second form is expressed in terms of the even (bonding) and odd (antibonding) 
operators 

Vje/oa — 

with e and o corresponding to + and — respectively. We will study the ground state rung- 
rung correlation function 

c(\k - j\) = {J k Jj) (4) 

Using DMRG techniques [|J we have calculated the rung-rung current correlation func- 
tion, Eq. (D, for a 2 x 32 ladder with J/t = 0.35 at fillings (n*) = 1 - x with x = 0.125 
and 0.19. The DMRG results are for a ladder with open-end boundary conditions and c(£) 
is obtained by averaging-over sites separated by £ = \k — j\. The calculations were done 
with the finite system DMRG method, keeping a maximum of 800 states per block, on a 
32 x 2 open system. The discarded weight was typically about 3 x 10~ 6 in the final sweep. 
The correlation measurements were checked by repeating some of the calculations keeping 
1000 states; the differences were quite small and would not affect any conclusions. In order 
to reduce the effect of the open boundaries, results for many different points with the same 
separation were averaged over to obtain each data point shown. As shown in Fig. 1(a) for 
x = 0.125, the rung-rung current correlations exhibit an incommensurate oscillation with a 
small amplitude and decay exponentially. Fig. 1(b) shows a semi- log plot of \c(£)\ versus i 
for J/t = 0.35 and x = 0.125 and Fig. 1(c) shows a similar plot for x = 0.19. The slope 
of the solid line in Fig. 1(b) is 1/2.8 corresponding to an exponential decay of c(£) with a 
correlation length £ = 2.8 lattice spacings. Similarly, as seen in Fig. 1(c), for x = 0.19 we 
find an exponential decay of the rung-rung current correlations with £ = 4 lattice spacings. 

Further insight into the nature of rung-rung current correlations on a 2-leg t-J ladder 
is provided within a bosonization framework. Here we pass to a continuum limit along the 
length of the 2-leg ladder and introduce the usual left- and right-moving fields 

<Mx) = e~ ik ^ LXa (x) + e lk ^ RXa (x) (5) 

with A = e and o corresponding to the even (bonding) and odd (antibonding) bands respec- 
tively with kp e and kp the fermi wave vectors for these bands. 
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Using this representation, the rung current operator, Eq. (Q), becomes 

J(X) = J (X) + ^Fe + k Fo )x (x) + h ^ (g) 



with 



J (x) = ft£[ e «<**.-**> (iPl a (x)^ Loa (x) - + h.c] (7) 



and 



J2k F = ttJ2(^Lea( X )^Roa(x) ~ lfj[ oa (x)^ Rea (x)^ . (8) 
a 

In the usual way, we represent the left and right moving fermion fields by left and right 
moving boson fields: 

IpL/RXa ~ e^ViUa (9) 

Then, introducing the dual canonical Bose fields, 

4>\a = 4>R\a + 4>L\a &Xa = 4>RXa ~ 4>L\a (10) 

we have, for example, for the first term in Eq. (§) 

^ Lea (x)tlj Roa (x) ~ e ^(~^ a +e ea +^ oa +e oa ) _ ^ 

In terms of the even and odd parity charge and spin Bose fields with A = e or o 

K = (0a T + M/V2 

ACT = (0AT-0Ai)/V2 (12) 
and their dual fields, we can write the up spin part of Eq. fllTD as 

VLt(z)VW^) ~ e^-*-,-*--^,-*^) . (13) 

Here, 

±p = (0 ep ± op ) /V2 (14) 



with similar relations for 4>± a , 9± P , and 9± a . Actually, this last transformation is not canon- 
ical when the even and odd bosons have different velocities. However, we will follow the 



standard practice |||6| of assuming that this velocity difference is irrelevant. In this way, 
the rung current operators, Eqs. (0) and (|8|) become 

J a (x) ~ z 5 1 e iw(kF e- kF ° )x e~ i ^- p+S2 *-°- Sie - p ~ Sl82e -° ] + h.c. (15) 

<Si,<5 2 =±l 

(where 5i labels the left or right term and 5 2 labels spin) and 

J2k F {%) ~ E 5 1 e t ^ [ ~ Sl ' t, -»~ 5lS2 ^ +e+ » +52e+ «] (16) 

5i,fo=±l 

(where now 5i labels the e — o or o — e terms and 52 labels spin) . 

Renormalization group analysis 0,0], based on the weak coupling 2-leg Hubbard model 
and earlier numerical results |7|-|9| on t-J ladders suggest that, for a realistic parameter range 
away from half-filling, the t-J model is in a "C1SO" phase in which both spin bosons are 
gapped and one of the charge bosons is gapped. When a boson is gapped, either ((f)) ^ or 
(6) 7^ 0. These expectation values must be specified to completely specify the phase which 
the system is in. The expected phase has 

(ftfar) + , (0_ p ) ± . (17) 

In this phase the uniform part of the pair correlation function has power-law decay and the 
2kp e + 2kp Q part of the density correlation function has power law decay (but not the 2kp e 
or 2k po parts), in agreement with the numerical results. 

We may replace factors like e l ^^~ p by their expectation values, so that: 

■hk F ~ Y, d^^ 1 ^-^^ . (18) 

Si,S 2 =± 

Correlation functions involving exponentials of decay exponentially since is conju- 
gate to the pinned dual field (Q- a )- Therefore, we expect that the term in the rung current 
correlation function which oscillates at 2kp will decay exponentially. Similarly, since J D de- 
pends exponentially on and 6L p and both of these are conjugate to pinned phases, these 
correlation function will also decay exponentially. 

Based upon this bosonization analysis, we expect that the asymptotic form of the rung- 
rung current correlations will oscillate with an incommensurate wave vector {kp + kp e ) = 
7r(l — x). Thus, for x = 0.125 the bosonization result gives 
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c{£) ~ e~ t/( cos(7tt£/8 + 0) (19) 

Using the DMRG results, £ = 2.8 for x = 0.125, Fig. 2 shows a comparison of Eq. ( |19|) with 
the DMRG data. 

We note that other phases have been suggested for the 2-leg Hubbard model and extended 
Hubbard models [|5l,|T0|| some of which could have power law decay for the 2k p component 
of the rung current correlation function. However, power law decay could only occur in a 
phase in which none of the boson fields dual to the ones appearing in Eq. (16) are pinned. 
This would mean that each boson field in Eq. (16) is either itself pinned or else is gapless. 
A possible C1S0 phase with this property has: 

(e +a )^o, (0- CT )^o, (<P- p )?o. (20) 

In fact, such a CI SO phase was predicted by Fabrizio in the Hubbard model for small t± 
and a small range of U near 7, based on a perturbative renormalization group analysis. Thus, 
it is possible that such a phase might also occur in the 2-leg t-J model or some generalization 
of it, for some range of parameters. However, Eq. (|20|) is not a sufficient condition for power 
law decay. Not only must these fields have non-zero ground state expectation values, but also 
these values must be such that a cancellation of the power-law part of the correlation function 
between the various terms in Eq. ( Jl6| ) does not occur. By explicitly writing pairing operators 
in the (f>/9, p/o", +/ — basis, used above, it can be checked that all pair correlations (singlet 
and triplet) have exponential decay in such a putative phase. We emphasize again that, based 
upon the d x 2_ y 2-Y\ke power law pairing correlations and density power law correlations which 
are observed for the t-J ladder for physically relevant J/t and doping regimes of interest, it 
would appear that the ground state of this model is characterized by pinned phases given 
by Eq. fllTD- In this ground state the current-current correlations decay exponentially. 
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FIG. 1. The current-current correlation function c(Z) on a 32 x 2 lattice calculated with DMRG, 
at dopings of x = 0.125 and x = 0.19, with J/t = 0.35. In (a), we show results on a linear scale 
for x = 0.125, and in (b), we show the same results on a semi-log scale. The solid line in (b) 
corresponds to a decay length of £ = 2.8. In (c), we show semi- log results for x = 0.19. The solid 
line corresponds to a decay length of £ = 4. 
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FIG. 2. (a) The current-current correlation function c(l) data from Fig. 1(a) is compared with 
the asymptotic form expected from bosonization. Here A = 0.0016, <fi = —2.0, and as in Fig. 1(b) 
the decay length is £ = 2.8 lattice spacings. 
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